Structural transitions in the 309-atom magic number Lennard-Jones cluster 
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The thermal behaviour of the 309-atom Lennard-Jones cluster, whose structure is a complete 
Mackay icosahedron, has been studied by parallel tempering Monte Carlo simulations. Surprisingly 
for a magic number cluster, the heat capacity shows a very pronounced peak before melting, which 
is attributed to several coincident structural transformation processes. The main transformation 
is somewhat akin to surface roughening, and involves a cooperative condensation of vacancies and 
adatoms that leads to the formation of pits and islands one or two layers thick on the Mackay 
icosahedron. The second transition in order of importance involves a whole scale transformation 
of the cluster structure, and leads to a diverse set of twinned structures that are assemblies of 
face-centred-cubic tetrahedra with 6 atoms along their edges, i.e., one atom more than the edges 
of the 20 tetrahedra that make up the 309-atom Mackay icosahedron. A surface reconstruction 
of the icosahedron from a Mackay to an anti-Mackay overlayer is also observed, but with a lower 
probability. 

PACS numbers: 61.46.Bc,36.40.Ei,36.40.Mr 
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I. INTRODUCTION 

Lcnnard- Jones (LJ) clusters have been widely stud- 
ied over the last few decades and have become a model 
system for understanding some of the structural, ther- 
modynamic and dynamic properties that are particu- 
lar to clusters. This choice is mainly due to the sim- 
plicity of the model, but it also reflects the possibility 
of making comparisons with experiments on rare gas 
clusters. For example, the ground-state structure of 
all LJ clusters with up to 1600 atoms has now been 
investigated! 1 : 2 : 3 : 4 : 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 : 13 It is well-established 
that, for N < 1000, LJ clusters follow an icosahedral 
pattern growth, showing magic numbers corresponding 
to complete Mackay icosahedra at TV = 13, 55, 147, 309, 
561, and 923. In between these magic numbers, the vast 
majority of the global minima correspond to a Mackay 
icosahedron core covered by an incomplete outer layer. 
The exceptions occur at or near to shell closings of al- 
ternative structural forms. Mostly, these correspond to 
Marks decahedra, but there is also one instance of a face- 
centred-cubic (fee) truncated octahedron (N — 38) 2,3 
and of a Leary tetrahedron (N = 98). 9 

There are two different ways in which a layer can 
grow on the surface of a Mackay icosahedron, which 
are termed Mackay and anti-Mackay growth modesi A 
Mackay icosahedron can be considered to be made up of 
twenty slightly distorted fee tetrahedra sharing a com- 
mon vertex. The Mackay growth mode continues the fee 
packing of the underlying tetrahedra, leading to the next 
Mackay layer, whilst the anti-Mackay growth mode in- 
volves sites that are in hexagonal close packing positions 
with respect to the underlying tetrahedra, thus intro- 
ducing a twin plane. The two overlayers are illustrated 
in Fig.n The growth of a new shell in LJ clusters usu- 
ally starts at the anti-Mackay sites, as this overlayer does 
not involve low-coordinate edge atoms. However, as the 
Mackay overlayer is associated with a higher surface den- 
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FIG. 1: Atomic positions for the two possible overlayers on 
a 147-atom Mackay icosahedron. 



sity, it has a lower strain energy and becomes energeti- 
cally favourable over the anti-Mackay above a given size. 

The thermal behaviour of LJ clusters has also been 
much ctiidicd . 14 ' 15 - 16,17,18,19,20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 - 28,29 - 30 It 
is known that many clusters do not retain the ground- 
state structure up to melting, but, instead, they un- 
dergo structural changes as the temperature is raised. 
For LJ clusters, solid-solid transitions prior to melting 
can be expected for those sizes with non icosahedral 
geometriesi 1 ?! 20 ' 21 : 25 : 26 : 27 : 28 : 30 Premelting effects associ- 
ated with structural changes at the surface that lead to 
enhanced surface diffusion are also commonii 6 ^ Of par- 
ticular interest to us here are the surface reconstructions 
which can occur for clusters with an incomplete Mackay 
overlayer and in which the outer layer changes from 
Mackay to anti-Macka y 24 : 29 This surface reconstruction 
is mainly driven by the greater vibrational entropy as- 
sociated with anti-Mackay structures. For a particular 
shell, the transition is first seen at the crossover size 
at which a Mackay structure becomes lowest in energy 
(e.g. at N = 31 for the second shell, and at N = 82 for 
the third shell), and also occurs at the subsequent sizes 
but at a temperature that increases with the cluster's 



2 



size. Therefore, there comes a size for which the melting 
temperature becomes lower than the temperature neces- 
sary for the surface reconstruction to occur and, so, the 
surface reconstruction only takes place in a limited size 
range, namely N = 31 — 3922i2li2S for the second Mackay 
layer and N = 82 - 140 29 for the third Mackay layer. 

As the growth of the fourth shell of the Mackay icosa- 
hedron follows the same pattern as the previous ones, it is 
very likely that the surface reconstruction again occurs. 
The results for the second and third layer seem to indicate 
that there might be a tendency for the size range over 
which the surface reconstruction takes place to increase 
with the number of shells. In particular, if the upper 
size for which this transition occurs continues to increase 
with the shell number, it might be that even the complete 
309-atom Mackay icosahedron could undergo such a sur- 
face reconstruction as well. This would be an unexpected 
result, since usually magic numbers clusters will show a 
single 'first-order-like' melting peak in the heat capacity, 
and 'premelting' and other transitions are more generally 
associated with incomplete shells. Indeed, Labastie and 
Whetten's classical paper on the size evolution of clus- 
ter melting towards the bulk first-order behaviour, was 
based on the Mackay icosahedral LJ clusters with 13, 55 
and 147 atoms^ However, there is some precedent for 
such a transition in a complete Mackay icosahedron. For 
Morse clusters, it has already been shown that such a 
surface reconstruction occurs prior to melting for N — 
561 and 923^2i although the parameter in the Morse po- 
tential that determines the range of the potential had 
a value for which anti-Mackay structures are stabilized 
relative to the LJ potential^ Therefore, a study of the 
melting behaviour of the 309-atom LJ cluster would serve 
both to test our hypothesis concerning a possible surface 
reconstruction and to better understand the thermody- 
namics of surface transitions in clusters. 



II. METHODS 

The potential energy of a cluster interacting with the 
LJ potential is given by 
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where e is the pair well depth and 2 1 / 6 ct is the equilibrium 
pair separation. 

The thermal behaviour of the LJ 309-atom com- 
plete Mackay icosahedron was studied using canonical 
Monte Carlo (MC) simulations with the parallel tem- 
pering exchange scheme (PT)iS The PT MC simula- 
tions consisted of 32 trajectories at temperatures rang- 
ing from 0.360 e/fc^ to 0.453 e/fcg. A preliminary study 
over a broader range of temperatures, from O.OlOe/fcs 
to 0.600 e/ks, showed that the region 0.360 e/ks - 
0.453 e/fce was wide enough to capture all the interest- 
ing features of the heat capacity curve. Due to the slow 



convergence, the system was equilibrated by performing 
7.2 xlO 7 MC cycles for each one of the trajectories, and 
the data to compute the heat capacity were collected 
over the following 3.06x10 s MC cycles. Similar conver- 
gence problems have already been encountered in previ- 
ous PT MC studies of the thermal behaviour of smaller 
LJ clustersi 20 i 21 i 26 i 27 i 28 i 29 Exchanges between two trajec- 
tories were attempted with a probability of 10%. All the 
temperatures were initialized with the Mackay icosahe- 
dron ground-state structure. Evaporation and fragmen- 
tation of the cluster was avoided by adding a hard wall 
of radius r + a, where r is the radius of the ground-state 
structure. The heat capacity and caloric curves were ob- 
tained from the simulation data using the multihistogram 
technique 



III. RESULTS 

In Fig. |21 we have plotted the canonical and micro- 
canonical caloric curves and the canonical heat capac- 
ity, as obtained from our simulations using the multihis- 
togram technique. There are two fairly sharp changes in 
the slope of the canonical caloric curve, which, in the mi- 
crocanonical ensemble, lead to two van der Waals loops, 
i.e. regions in which the caloric curve exhibits a negative 
slope and the heat capacity is negative. Coincident with 
these features in the caloric curves, at the temperatures 
0.390 e/ks and 0.417 e/ks, the canonical heat capacity 
curve shows two clearly distinct peaks. The height of the 
first peak is roughly half the height of the second one, but 
the former is slightly broader. The latent heat of these 
transitions was computed as the area under each peak in 
the heat capacity curve. It was found that the latent heat 
associated with the first peak is approximately 0.67 l m , 
where l m = 0.29/iVe is the latent heat of the second peak, 
which, as we will show below, signals the melting of the 
cluster. These results indicate that the 309-atom LJ clus- 
ter undergoes a major structural transition before melt- 
ing. The unusual presence of two van der Waals loops in 
the microcanocial caloric curve, as well as the two strong 
peaks in the canonical heat capacity, suggest that this 
transition might be, as with the cluster melting transi- 
tion, the finite-size analogue of a first-order transition. 15 
In the regions where the microcanonical heat capacity is 
negative, the canonical probability distribution for the 
energy is also bimodali^i 

The structural transformations associated with the two 
peaks in the heat capacity curve were analyzed by per- 
forming systematic quenches starting from configurations 
generated in each of the 32 trajectories employed in 
our calculations^ 5 - Two useful properties to character- 
ize the minima are their potential energy and geomet- 
ric mean vibrational frequency v. The main contribu- 
tions to v come from atoms which are nearest neighbours. 
For a pair potential, such as LJ, the value of v reflects 
both the number of nearest neighbours, and the aver- 
age nearest-neighbour distance, the latter because the 
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FIG. 2: (a) Canonical (solid line) and microcanonical (dashed 
line) caloric curves and (b) canonical heat capacity for the 
309-atom LJ cluster. 



second derivative of the LJ potential decreases with dis- 
tance. Strained structures, which inevitably have a vari- 
ety of nearest-neighbour distances, have a larger average 
nearest-neighbour distance, because of the anharmonic 
nature of the LJ well, and hence have smaller values of 
v. 

Plotting v against the potential energy for the minima 
we obtained shows that they can be naturally divided 
into a number of types, which, as we will see, reflect their 
different structural features (see Fig.|3£a)). Specifically, 
the diagram can be divided into four well-defined zones, 
with only a relatively small number of minima near to 
the boundaries. A closer inspection of zone 2 shows that 
this zone can be further divided into two subzones (2A 
and 2B) on the basis of their frequencies, but the division 
in this case is less clear cut. 

To determine the structural types associated with 
the different zones in Fig.[3Ja) we visually inspected 
a large number of minima. To facilitate this task, 
wc also performed common neighbour analyses of the 
structuresj2£i2L2& and so we were able to identify atoms 
that had a particular local coordination environment. It 
was particularly noteworthy that all the minima in zones 
1, 2 and 3 are ordered in the sense that all the atoms in 
the interior of these clusters have twelve nearest neigh- 
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FIG. 3: (a) (Colour online) A scatter plot of the geometric 
mean vibrational frequency against the potential energy for 
the minima collected from the quenches at all temperatures. 
Five distinct zones can be distinguished in this diagram, and 
the boundaries that we have chosen to delineate these zones 
are explicitly shown in the plot. The unit of frequency is 
(e/mcr 2 )' 1,/2 - > . (b) The percentage of quenches that lead to 
minima in each of the five zones as a function of temperature. 



bours, whose local coordination shell can be character- 
ized as icosahedral, decahedral, fee or hep^ Typical iso- 
mers representative of each zone are depicted in Fig.^J 

The lowest energy region (Energy < — 1972 e), desig- 
nated as zone 1, corresponds to the global minimum, 
the complete Mackay icosahedron, and structures based 
upon it, in which there are some vacancies at the 
low-coordinate vertex and edge sites, and the resulting 
adatoms are dispersed across the surface of the cluster, 
perhaps forming adatom clusters of a few atoms (see 
structure (1) in Fig. As the number of vacancies 

increases, the potential energy of the minima increases, 
and their mean vibrational frequency decreases due to 
the lowered number of nearest neighbours. 

At energies in the interval — 1972 e to — 1915 e, there 
are two clearly separate regions differentiated by their 
vibrational frequency. As mentioned above, the higher 
frequency region, zone 2, is further divided into two sub- 
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zones designated as 2 A and 2B. The isomers found in 
zone 2A are Mackay icosahedra in which a pit is formed 
on one side of the cluster, whilst an adatom island is 
formed on the opposite side (see structures (2A) I and II 
in Fig.0J). The islands prefer to occupy anti-Mackay sites, 
and both the pits and islands can be one or two layers 
thick. These structures are therefore similar to those of 
zone 1, but have a significantly larger density of adatoms 
and vacancies. Consistent with this, the zone 2A minima 
form a downwardly sloping band in Fig.^a), that can 
be seen to be an extension of the band due to the zone 
1 minima. However, there is a clear gap between these 
two bands at the boundaries between zones 1 and 2A, 
because minima with an intermediate number of surface 
defects are less likely to be sampled. 

It used to be common to refer to a Mackay icosahe- 
dron as a multiply-twinned particle^ because they are 
made of twenty somewhat strained fee tetrahedra shar- 
ing a common vertex. The twin planes correspond to the 
internal faces of these tetrahedra and which are common 
to two tetrahedra. For the 309-atom Mackay icosahedron 
there are five atoms along the edges of the fee tetrahedra. 
Other multiply twinned structures that are also made up 
of aggregates of fee tetrahedra are possible, such as the 
Marks decahedral and the Leary tetrahedron. 9 

After an extensive investigation of the structures of 
the minima in zone 2B, in which the common neighbour 
analysis was particularly revealing, we found that the 
common feature of this diverse set of structures was that 
they are aggregates of fee tetrahedra with 6 atoms per 
edge. As the fee tetrahedra are larger than for structures 
based on the 309-atom Mackay icosahedron, the struc- 
tures are less strained, and so have larger vibrational 
frequencies. Three different examples of the possible ar- 
rangements of the fee tetrahedra are shown in Fig.^J For 
structure (2B) I there is a single complete fee tetrahedron 
at the centre of the cluster and it is related to a Leary 
tetrahedron £ Structure (2B) II has two twinned fee tetra- 
hedra at its centre. Finally, structure (2B) III has five fee 
tetrahedra (some incomplete) sharing a common edge, 
and so it has a decahedral structure at its core. In all 
these structures, these central tetrahedra are covered by 
overlayers that begin the growth of new tetrahedra that 
are twinned with respect to the core. 

The lowest vibrational frequency region at intermedi- 
ate energies, zone 3, is occupied by isomers with a 147- 
atom Mackay icosahedral core, but where the outer layer 
has undergone a surface reconstruction from a Mackay to 
an anti-Mackay overlayer. Because of the greater strain 
associated with these structures, they have a lower mean 
vibrational frequency. As mentioned in the introduc- 
tion, a similar surface reconstruction has been observed 
for some smaller LJ clusters with incomplete Mackay 
layers^S anc ; f or larger complete Mackay icosahedra in- 
teracting with a long-ranged Morse potential^ As seen 
in Fig.|3|(a), this surface reconstruction is partly driven 
by the higher vibrational entropy of the anti-Mackay 
overlayer. In addition, covering the 147-atom Mackay 



icosahedron with a complete anti-Mackay overlayer gives 
rise to a 279-atom structure, and so there is also a con- 
figurational entropy gain associated with the possible ar- 
rangements of the 30 extra adatoms. 

Note that the edges of the underlying fee tetrahe- 
dra around which the above anti-Mackay structures are 
formed have only four atoms along their edges. There- 
fore, we can conclude that the three regions at interme- 
diate energies correspond to structures built around as- 
semblies of fee tetrahedra, but, depending on the length 
of the edges of these tetrahedra, and hence the associated 
strain, the isomers fall in different regions of frequency. 

Finally, the high energy region, zone 4, is associated 
with completely melted structures, in which both the core 
and surface atoms are disordered. There is a variety of 
coordination numbers associated with the interior atoms 
in these structures. Furthermore, the common neighbour 
analysis shows that the vast majority of the atoms do not 
have an ordered local coordination shell. Therefore, zone 
4 is assigned to the liquid state. 

An analysis of the quench results is presented in 
Fig-El(b) showing the percentage of quenches that lead 
to minima in the different zones. At temperatures below 
that for the first peak almost all isomers found corre- 
spond to Mackay icosahedra in which a number of vacan- 
cies and adatoms are formed on the cluster's surface. The 
vacancies mainly appear at the low-coordinate vertex and 
edge sites, whilst the adatoms diffuse around the surface, 
as seen previously for smaller Mackay icosahedra>iSili At 
temperatures coincident with the first peak there is a 
fairly sharp transition from the Mackay icosahedron with 
a small number of defects to structures which have under- 
gone some kind of structural transition. Surprisingly, the 
Mackay to anti-Mackay reconstruction represents only a 
minor contribution to the heat capacity, with only up to 
around a 15-20% probability of finding the cluster with 
an anti-Mackay reconstructed surface. Instead, two other 
transitions dominate in this range of temperatures. 

In the main transition, which represents a 50-55% of 
the total transition, adjacent atoms from the outer shells 
of the Mackay icosahedron vacate their sites, leaving 
behind pits one or two layers thick, and the resulting 
adatoms form islands above the still intact parts of the 
fourth shell. As mentioned above, below the transition, 
the vacancies and adatoms associated with the Mackay 
icosahedron are dispersed across the surface and their 
number gradually increases with temperature. However, 
as the surface density of these defects increases there 
comes a point where there is a jump in the number of 
vacancies and adatoms, and they cluster together to form 
islands and pits. This process is driven by the (effective) 
attraction between the adatoms and vacancies. Once 
a number of vacancies are formed at vertex and edge 
sites, the neighbouring atoms also lower their coordina- 
tion number and, therefore, it becomes more favourable 
to remove more vertex and face atoms in positions adja- 
cent to an already existing vacancy or pit. This transition 
is somewhat akin to the pre-roughening and roughening 



5 




(2A)II 




E=-2007.2i90 
v =2.256 



E=-i997.7675 
v =2.239 



(2B)I 



E=-i966.3657 
v =2.186 



(2B) II 



E=-i964.7I86 
v =2.190 





(2B)III 



E=-1965.7417 
S =2.213 



E=-i959.0320 
1) =2.196 




E=- 1963.41 05 
1) =2.210 



E=-l 956.5963 
l. 2.116 



FIG. 4: (Colour online) Structure of the global minimum and of some of the isomers representative of the zones shown in Fig. 
3(a). The energy and frequency are given in reduced units. To illustrate how the structures can be viewed as assemblies of fee 
tetrahedra, to the right of some of the complete structures, we have also depicted the network formed by those atoms that lie 
along the edges of these tetrahedra. 



transtionS)** both of which are known to occur for the 
bulk LJ {111} surfaces4M4 

The second transition competing with the Mackay to 
anti-Mackay surface reconstruction, is not a surface tran- 
sition, but involves a transformation of the whole struc- 
ture of the cluster. The resulting structures are multiply- 
twinned particles that are assemblies of fee tetrahedra 
with six atoms along each edge, one atom more than 
the tetrahedra that make up the 309-atom Mackay icosa- 
hedron. This transition must be driven by an increase 
in the configurational entropy, as these structures have 
lower vibrational entropies. Indeed, this conformational 
entropy is reflected in the diversity of structures that we 
observed (Fig.^J. The probability that these structures 
are adopted reaches a maximum of about 30% midway 
between the two heat capacity peaks. 

Finally, as the second heat capacity peak is approached 
and the cluster begins to melt, this array of different 
structures becomes less probable in favour of disordered 
structures typical of the liquid state. 



IV. CONCLUSIONS 

The existence of a well-defined structural transition in 
the 309-atom LJ cluster prior to melting suggests that, 
as we anticipated in the introduction of this paper, there 
is a tendency for the size range for which LJ clusters 
with Mackay overlayers undergo structural transforma- 
tions prior to melting to increase as the number of shells 
in the cluster increases. For the fourth icosahedral shell, 
we expect these transitions to begin at around N = 169, 
the size at which a Mackay overlayer becomes energeti- 
cally more favourable than the anti-Mackay, and to con- 
tinue up to, and probably beyond, the next complete 
Mackay icosahedron at TV = 309. This compares to 
N = 31 - 39 and N = 82 - 140 for the second and third 
icosahedral shells. However, the unexpected aspect of our 
results was the nature of the transition that occurred and 
the size of the heat capacity peak that resulted. Contrary 
to our predictions, the Mackay to anti-Mackay surface re- 
construction was only a minor contribution to the heat 
capacity peak. Instead, two other transitions had a more 
significant effect. In particular, the dominant contribu- 
tion is from a surface transition involving the condensa- 
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tion of vacancies and adatoms that leads to the formation 
of pits and islands one or two layers thick on the Mackay 
icosahedron. Also fairly important was another transi- 
tion that leads to a complete loss of the Mackay icosa- 
hedral structure, and the adoption of multiply-twinned 
structures based upon packings of fee tetrahedra with six 
atoms along each edge. By contrast, the twenty strained 
fee tetrahedra that make up the 309-atom Mackay icosa- 
hedron have five atoms along each edge. 

The substantial latent heat, the van der Waals loop in 
the microcanonical caloric curve, and the bimodal dis- 
tribution for the energy in the canonical ensemble, all 
suggest that the transition at T=0.390e//cB is the finite 
size analogue of a first order phase transition. But does 
it have any analogues in bulk systems? The two more mi- 
nor transitions that contribute to the peak clearly have 
no bulk analogue, as they are specific to Mackay icosa- 
hedra and multiply-twinned particles, both of which will 
not be energetically competitive in the bulk limit. How- 
ever, as we have already mentioned, the cooperative con- 
densation of the vacancies and adatoms to form pits and 
islands, respectively, has some similarities to the bulk 
pre-roughening and roughening transitions 42 of the LJ 
fee {111} surfaces .42^4 All involve the loss of the uni- 
form, flat character of the perfect {111} surfaces through 
the formation of islands and pits, whilst retaining the ba- 
sic lattice structure. However, the pre-roughening tran- 
sition leads to a disordered flat surface, where the height 
fluctuations are restricted to be of one layer thickness, 
whereas, for the clusters, islands and pits are seen that 
can be two layers high or deep. The surface roughen- 
ing transition is defined as the temperature at which a 
height- height correlation function diverges with distance. 
However, such a definition cannot be straightforwardly 
applied to our cluster, because the {111} facets are so 
small. Furthermore, the cooperative nature of the tran- 
sition is also related to the small size of the facets, since 
it allows the removal of whole faces, nucleated by the 
preferential generation of vacancies at the vertices and 



edges. Therefore, it would premature to say that the 
transition in our clusters is the finite size analogue of ei- 
ther the pre-roughening or roughening transitions, and 
it would be interesting to see how the transition evolves 
with cluster size. 

Our results have also shown that even for a very strong 
magic number, such as the 309-atom LJ cluster, the heat 
capacity curve can exhibit strong features prior to melt- 
ing. This is an important finding, because usually it is 
expected that magic numbers will show only one well- 
defined melting peak in the heat capacity curve, whilst 
premelting effects and other structural transitions are 
more generally associated with clusters with incomplete 
shells. For example, following the scaling arguments of 
Labastie and Whetten, which were based on the first 
three LJ Mackay icosahedra^ it would have been ex- 
pected that the 309-atom cluster would exhibit a single 
heat capacity peak that was narrower and higher than 
for the 147-atom cluster, hence smoothly converging to- 
wards the bulk limit. Our results have shown that this 
an oversimplified view of the problem and that many sur- 
prising features can appear along the path towards the 
bulk limit. Indeed, one should expect structural transi- 
tions prior to melting to be the norm for larger LJ clusters 
up to very large sizes. Firstly, larger Mackay icosahedra 
would be expected to show similar transitions to those 
found here. Secondly, even at sizes when decahedral 
and then fee structures become lowest in energy, there 
are likely to be solid-solid transitions from decahedral to 
icosahedral and at larger sizes from fee to decahedral, 
driven by differences in vibrational entropy^ 
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